HYPERBOLIC RELAXATION OF REACTION DIFFUSION EQUATIONS WITH 
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Abstract. Under consideration is the hyperbolic relaxation of a semilinear reaction-diffusion equation, 

euu + ut - Au + f(u) = 0, 
on a bounded domain SlcS 3 , with e £ (0, f] and the prescribed dynamic condition, 

d n u + u + ut = 0, 

on the boundary Y := dQ. We also consider the limit parabolic problem (e = 0) with the same 
dynamic boundary condition. Each problem is well-posed in a suitable phase space where the global weak 
solutions generate a Lipschitz continuous semiflow which admits a bounded absorbing set. Because of 
the nature of the boundary condition, fractional powers of the Laplace operator are not well-defined. The 
prccompactness property required by the hyperbolic semiflows for the existence of the global attractors 
is gained through the approach of 1441 . In this case, the optimal regularity for the global attractors 
is also readily established. In the parabolic setting, the regularity of the global attractor is necessary 
for the semicontinuity result. After fitting both problems into a common framework, a proof of the 
upper-semicontinuity of the family of global attractors is given at e = 0. Finally, we also establish the 
existence of a family of exponential attractors. 
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1. Introduction 

Let fl be a bounded domain in M 3 with boundary T := d£l of class C 2 . We consider the hyperbolic 
relaxation of a semilinear reaction diffusion equation, 

£u tt + ut-Au + f(u) = 0, (1.1) 

in (0, oo) x f2, where e € [0, 1]. The equation is endowed with the dynamic boundary condition, 

d n u + u + u t =0, (1.2) 

on (0, oo) x r, and with the initial conditions, 

m(0, x) = uq(x), ut(0, x) = Ui(x) in Cl. (1-3) 
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For the nonlinear term /, we assume, / € C 2 (R.) and that there is a constant, I > 0, such that, for all 
s£l, the following growth and sign conditions are satisfied, 

l/"0OI<*(i + M), (i.4) 

and 

f(s) 

liminf — > -A, (1.5) 

|s|->oo S 

where A > is the best Sobolev/Poincare type constant 



A / u'dx< / \Vu\ J dx+ / u'dS. (1.6) 
Jn Jn Jr 

Finally, assume that there is $ > such that for all s € M., 

f'(s) > -0. (1.7) 

Notice that the derivative / = F of the double- well potential, F(u) = jU 4 — ku 2 , k > 0, satisfies 
assumptions (|1.4[) , (|1.5|) and (|1.7[) . The first two assumptions made here on the nonlinear term, (|1.4I) 
and (|1.5p . are the same assumptions made on the nonlinear term in |16j . |41j and |51j . for example 
([H] additionally assumes /(0) = 0). The third assumption (fT77|) appears in [13]. [2B"] . [27] and 
the bound is utilized to obtain the precompactness property for the semiflow associated with evolution 
equations when dynamic boundary conditions present a difficulty (e.g., here, fractional powers of the 
Laplace operator subject to ([1.2)1 arc undefined). It is worth mentioning that ([1.5)1 can be also replaced 
by a less general (but still widely used in the literature) condition 

liminf / (s) > — A 

in which case, (|1.7)) is automatically satisfied. Furthermore, assumption (|1 .4[) implies that the growth 
condition for / is the critical case since Q C R 3 . Such assumptions are common when one is investigating 
the existence of a global attractor or the existence of an exponential attractor for a partial differential 
equation of evolution. 

Of course, when (|l.ip is equipped with Dirichlet, Neumann or periodic boundary conditions, (]1.6I) 
simplifies. Moreover, if (|1.1[) is equipped with a Robin boundary condition, then an estimate like ([1.6) 1 
holds, but A possesses an explicit description as the first eigenvalue of the Laplacian with respect to 
the Robin boundary condition. The relation between the dynamic condition ([1.2)1 with the acoustic 
boundary condition is discussed below. The hyperbolic equation (|1.1[) is a well-known nonlinear wave 
equation motivated from (rclativistic) quantum mechanics (cf., e.g [3"l ll5[[3"6ll49j ). However, as mentioned, 
most sources study the asymptotic behavior of ([1.11 1 with a static boundary condition such as Dirichlet, 
Neumann, periodic or Robin. One of the goals of this paper is to extend some results concerning the 
asymptotic behavior of ([1.1)1 , now with the dynamic boundary condition ([1.2)1 . The corresponding linear 
case for (|1.1)) - (|1.3|) is treated in [46]. The existence of the global attractor for a linear damped wave 
equation with a nonlinear dynamic boundary condition is considered in [53] . More general systems, with 
supercritical nonlinear sources on both the interior and the boundary, are considered in [3J [S] QUI US H2] • 
These contributions mainly devote their attention to issues like, Hadamard local wellposedness, global 
existence, blow-up and non-existence theorems, as well as estimates on the uniform energy dissipation 
rates for the appropriate classes of solutions. We also refer the reader to [13] for a unified overview of 
these results. 

Our main goal is to compare the hyperbolic relaxation problem (|1.1[) - (I1.3I) with that of the limit 
parabolic equation where, for e — 0, we have the reaction-diffusion equation, 

ut-Au + /(«) = (), (1.8) 

in (0,oo) x 17, with the dynamic boundary condition, 

d n u + u + u t =0, (1.9) 

on (0,oo) x r, and the initial conditions 

u(0,x) = uq(x) in Q, u(0,x) = Jq(x) on T. (1-10) 
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For the sake of simplicity, we shall restrict our attention only to linear boundary conditions of the form 
(|1.9j) even though our framework can easily allow for a complete treatment of nonlinear dynamic boundary 
conditions (see Remark GuU cf. also [14], [23], [27]). 

Because of its importance in the physical sciences and the development of mathematical physics, the 
reaction-diffusion equation (|1.8p and its asymptotic behavior are well-known to the literature. Many of 
the books referenced above contain a treatment on the parabolic semilinear reaction-diffusion equation 
(|1.8[) with static boundary conditions. In particular, the Chaffcc-Infante reaction-diffusion equation, with 
f(u) = u 3 — ku, k > and Dirichlet boundary conditions can be found in [57J Section 11.5]. A discussion 
on the structure of the associated global attractor can also be found there. Additionally, the Chaffce- 
Infante equation and its hyperbolic relaxation, again with Dirichlet boundary conditions, are discussed 
in H3 Chapters 3-5]. 

Recently there has been a great amount of research taking place in the area of partial differential 
equations of evolution type, subject to dynamic boundary conditions. Boundary conditions of the form 
(|1.9p arise for many known equations of mathematical physics. This can especially be seen by the many 
applications given to heat control problems, phase-transition phenomena, Stefan problems, some models 
in climatology, and many others. Without being too exhaustive we refer the reader to |2U d5J HH] for 
more details about the system (|1.8[) - (|1.10l) and a more complete list of references. A version of equation 
(|1.2j) . but with nonlinear dissipation on the boundary, already appears in the literature, we refer to 
[lH [TTl [18] . There, the authors are able to show the existence of a global attractor without the presence 
of the weak interior damping term u t , by assuming that / is subcritical. One motivation for considering 
a boundary condition like (|1.2p comes from mechanical considerations: there is frictional damping on 
the boundary T that is linearly proportional to the velocity u t . In [51j . the convergence, as time goes 
to infinity, of unique global strong solutions of (|1.1[) - (|1.3[) to a single equilibrium is established provided 
that / is also real analytic. Note that the set of equilibria for (jl.l[l - (jl.3[l may form a continuum so that, 
in general, guaranteeing this convergence is a highly nontrivial matter. The second motivation comes 
from thermodynamics. Suppose that we want to consider heat flow in a metal. The standard derivation 
of the heat equation is always based on the idea, that "heat in equals heat out" over a region Q. But 
the classical approach ignores the contribution of heat sources located on the boundary T. by taking into 
account only heat sources/sinks which are present inside the region (in our case, —f(u) is treated as a 
source within fi). A new derivation of the heat equation in the presence of heat sources/sinks located 
at r, assuming the Fourier law of cooling states (i.e., the heat flux is directly proportional to the 
temperature gradient, q = — Vu) was given in [31] . and it has lead to the precise formulation of the 
system in (jl .8[) - (|l . 10[) . However, the derivation in [31] suffers from an important drawback which cannot 
be ignored: initial perturbations in (|1.8[) propagate with infinite speed. This means that the presence of a 
heat source located at T is instantaneously felt by all observers in f2, no matter how far away from T they 
happen to be. This behavior can be traced to the "parabolic" character of Fourier's law. Thus, in many 
relevant phenomena the system (|1.8I) - (|1.10[) can become a bad approximation (sec, e.g., [T], [35], for many 
examples). In order to overcome these problems, a generalization of the standard Fourier law must be 
considered, leading to a new formulation for which the heat flux q obeys the so-called Maxwell-Cattaneo 
heat conduction law: 

ed t ~t + t = -Vu, (1.11) 

in (0, oo) x CI. Note that the Fourier law is obtained from (|1.11[) when e = 0. This expression for the heat 
flux ~<f leads to the hyperbolic equation (|1.1[) . which entails that u propagates at finite speed. It is also 
worth mentioning that one can write p. lip in the equivalent form of 

f°° 1 
lt(t,x) = - 9 e (t- s) Vu(s,x)ds, 9 e (t) := -e~i . (1.12) 
Jo £ 

This points to a situation in which the (past) thermal memory of the material plays a role, but its 
relevance goes down quickly as we move to the past. Finally, it may be worth mentioning that the form 
of flux if assumed in (|1.12p . in which e is assumed to be a generic memory kernel, also yields the 
following problem: 

/>oo 

u t = / e e (t-s)(Au(s)-f(u(s)))ds. (1.13) 
Jo 
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In this case, O e (s) = e~ 1 Q(s/e) and 9 : (0, oo) — > (0, oo) is a given (smooth) summable and convex 
(hence decreasing) relaxation kernel. A complete treatment of equation (|1.13[) , endowed with the dynamic 
boundary condition (|1.9j) . will be the subject of further investigation in the future. 

It may also be interesting to note that the dynamic boundary condition given in (|1.2[) can be recovered, 
in some sense, from the linear acoustic boundary condition, 

J m5 u + 5 t + 5 = -u t on (0,T)xT , . 

\ d n u = 6 t on (0,T) x T. 1 ' 

Here the unknown S = 8(t, x) represents the inward "displacement" of the boundary F reacting to a 
pressure described by ~u t . The first equation (|1.14lh describes the spring-like effect in which T (and S) 
interacts with — u t , and the second equation (|1.14[) -? is the continuity condition: velocity of the boundary 
displacement S agrees with the normal derivative of u. Together, (|1.14j) describes T as a locally reactive 
surface. The term m = m(x) represents mass, so in a massZess system, the inertial term disappears. In the 
case when 5 can be modelled by u near the boundary; i.e., if <5 ~ u near T, then we arrive at the boundary 
condition described by (|1.2p . In applications, the unknown u may be taken as a velocity potential of 
some fluid or gas in VL that was disturbed from its equilibrium. The acoustic boundary condition was 
rigorously described by Beale and Rosencrans in [6] and [7]. Various recent sources investigate the wave 
equation equipped with acoustic boundary conditions, [19j [26] 021 [50] . However, more recently, it has 
been introduced as a dynamic boundary condition for problems that study the asymptotic behavior of 
weakly damped wave equations, see [23] and |48j . 

The aim of this paper is to extend the asymptotic results for dissipative wave equations (|1.1[) and 
reaction-diffusion equations (|1.8[) with the dynamic boundary condition (|1.2[) , in terms of a perturbation 
problem, and ultimately discuss the continuity of the attracting sets generated by these problems. Due 
to the nature of the boundary condition imposed for the model problem , we are unable to prove the 
existence of global attractors for the hyperbolic relaxation problem through the compactness argument 
which is typical for damped wave equations with static boundary conditions, such as Dirichlet, Neumann, 
periodic, or Robin boundary conditions (cf. e.g. [40j [49j [52| ) . The problem arises from our lack to define 
fractional powers of the Laplacian with respect to the boundary condition (|1.2[) . This situation takes 
place because of the permanence of the Ut term on F, which in turn means the "Laplacian" is not self- 
adjoint. Thus, for example, the model problem does not enjoy an explicit Poincare inequality found with 
a Fourier scries, nor the existence of a local weak solution found with a typical Galcrkin basis. Local 
solutions will be sought with semigroup methods that rely on monotone operators techniques as in |16j . 
Then estimates are applied to extend the local solutions to global ones and the existence of an absorbing 
set is determined. For the hyperbolic relaxation problem ([l.l[) - (|1.3jh we obtain the relatively compact 
part in the decomposition of the solution by following the approach in [44| . 

The main novelties of the present paper with respect to previous results on the damped wave equation 
(|I.ip are the following: 

• We extend the results on the existence of global attractors {^ e } eG ^ ^ for the damped wave equa- 
tion with a critical nonlinearity and a "dynamic" boundary condition instead of the usual 
Dirichlet boundary condition (see, e.g., [33], [34]). This is achieved through the decomposition 
method exploited in [44] which allows us to establish that A £ has also optimal regularity (see 
Theorem Ell]). 

• We show that a certain family \ A e \ of compact sets, which is topologically conjugated to 

I J eS[0,l] 

{.A e } e6 [ ji in a precise way, is also upper-semicontinuous as e goes to zero. Roughly speaking, we 

show that these sets A e converge to the "lifted" global attractor Ao associated with the parabolic 
problem. The argument utilizes the sequential characterization of the global attractor (cf., e.g. 
[40] Proposition 2.15]). The main difficulty comes from the fact that the phase spaces for the per- 
turbed and unperturbed equations are not the same; indeed, solutions of the hyperbolic problem 
are defined for (uq, u\) E H s+1 (O) x H s (H), s <E {0, 1}, while solutions of the parabolic problem 
make sense only in spaces like L 2 (fi) x L 2 (T) and H s+1 (fl) x H s+1 ^ 2 (T) , respectively (see 
(|1.10p ). Thus, previous constructions obtained for parabolic equations with Dirichlet boundary 
conditions cannot be applied and have to be adapted. 
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• We prove the existence of a family of exponential attractors {.M e } , e € (0, 1], which entails that 
A e is also finite dimensional even in the critical case. We recall that the same result was shown 
in |16| for the wave equation (i.e., (|1.1|) without any damping in Q) subject to the boundary 
condition (|1.9j) only in the subcritical case. Unfortunately, we are unable to show that this 
dimension is uniform with respect to e > as s goes to zero. Some other open questions are 
formulated at the end of the article. 

The article is organized as follows. The limit (e = 0) reaction-diffusion problem is discussed in Section 
[2J The section is mostly devoted to citing the already known main results of the parabolic problem: 
the existence and uniqueness of global solutions in an appropriate phase space (see Theorem I2.3|) , the 
definition of the (Lipschitz) semiflow, the existence and regularity of the global attractor (see Theorem 
12. 6p . Section [3] contains our treatment of the hyperbolic relaxation problem, for all e £ (0,1]. We discuss 
the existence and uniqueness of solutions defined for all positive times in Section [3.21 (see Theorem 13. 6j) . 
The solutions generate a semiflow on the phase space, and thanks to the continuous dependence estimate, 
we know that the semiflow is locally Lipschitz continuous. The existence of a bounded absorbing set is 
also shown (see Lemma f3. 10|) . The global attractor and its properties are established in Scction [3~3l while 
the upper-semicontinuous result is established in Section [33] The existence of exponential attractors for 
the hyperbolic problem is presented in Section [4] The statement of a Gronwall-type inequality, used 
frequently in the estimates, is included in the Appendix. 



In this short section, we recall some results for the limit parabolic problem (|1.8l) - (|1.10p . i.e., (|1.1I) - (|1.3I) 
with e = 0. Unlike the hyperbolic problem, a full general treatment of the limit parabolic problem with 
dynamic boundary conditions already appears in the literature (cf., e.g., [24] [25] [29] [39] and references 
there in); in particular, this section will summarize some of the main results from [25j . It should be 
noted for the interest of the reader that all formal calculations made with the weak solutions of the 
parabolic problem can be rigorously justified using the Galerkin discretization scheme that appears, for 
instance, in [55] Theorem 2.6]. Indeed, it is through the use of the Galerkin approximations that the 
existence of weak solutions for the parabolic problem is shown. The solution operator associated with 
the parabolic problem generates a locally Lipschitz continuous semiflow on the appropriate phase space. 
We also know that this semiflow admits a connected global attractor that is bounded in a more regular 
phase space. It follows that solutions, when restricted to the global attractor, are in fact strong solutions, 
exhibiting further regularity that will become essential when we later consider the continuity properties 
of the family of global attractors produced by the hyperbolic relaxation problem (e > 0) and the limit 
parabolic problem [e = 0). 

We need to introduce some notations and definitions. From now on, we denote by || ■ |, || ■ \\k, the 
norms in L 2 (Vl) 1 H k {Q), respectively. We use the notation (■,■) and (v)fc to denote the products on 
L 2 (Q) and H k (il), respectively. For the boundary terms, || • \\l 2 (t) an d (•, -)_L 2 (r) denote the norm and, 
respectively, product on L 2 (T). We will require the norm in H k (T), to be denoted by j| • ||ijfc(r), where 
k > 1. The L p (£l) norm, p £ (0, oo], is denoted | • \ p . The dual pairing between i? 1 (fi) and its dual 
(i/ 1 (il))* is denoted by (u,v). We denote the measure of the domain f2 by |f2|. In many calculations, 
functional notation indicating dependence on the variable t is dropped; for example, we will write u in 
place of u(t). Throughout the paper, C > will denote a generic constant, while Q : R+ — > R+ will 
denote a generic increasing function. All these quantities, unless explicitly stated, are independent of e. 
Further dependencies of these quantities will be specified on occurrence. 

The following inequalities are straight forward consequences of the Poincare type inequality (|1.6[) and 
assumptions (|1.5|) and (|1.7p . From (|1.5|) it follows that, for some constants fi £ (0, A] and c\ = ci(/, |f2|) > 
0, and for all £ £ H 1 ^), 



2. The limit parabolic problem 



</(0,0>-(a-m)II£II 2 -ci 



(2.1) 
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Let F(s) = J„ s f(a)da. For some constant C2 = C2(/, \ > 0, and for all £ G 

F(0d^>-^|ia 2 -c 2 

(2.2) 

See [TSJ page 1913] for an explicit proof of (|2.2[) . The proof of (|2.1j) is similar. Finally, using (|1.7[) and 
integration by parts on F(s) = J f(a)da, we have the upper-bound 

F(0dx<(f(t),0 + iu\\ 2 



<(f(0,0 + ^U\\l 



(2.3) 



The natural energy phase space for the limit parabolic problem (|1.8j) - (|1.10[) is the space 

Y = L 2 (fl) x L 2 (r), 

which is Hilbert when equipped with the norm whose square is given by, for all C = (u, 7) G Y, 

\\a 2 Y--=\\u\\ 2 + h\\h ir) . 

It is well-known that the Dirichlet trace map tr^ : C°° (QJ — ^ C°° (r) , defined by tro (w) = U|r extends to 
a linear continuous operator tin '■ H r (O) — ^ fj r - 1 / 2 (r) , for all r > 1/2, which is onto for 1/2 < r < 3/2. 
This map also possesses a bounded right inverse tr^ 1 : H r ~ 1 l 2 (r) — > i/ r (f2) such that tro (tr^, 1 -)/)) = x/>, 
for any tp G H r - 1 / 2 (r). Identifying each function -0 € C (0) with the vector V = (t/>, trpj (i/>)) € 
C (fi) x C (r), it follows that C (fi) is a dense subspace of Y = L 2 (fl) x L 2 (r) (see, e.g., [43l Lemma 
2.1]). Also, we introduce the subspaces of H r (f2) x iJ r_1 / 2 (r), for every r > 1/2, 

V := {(u >7 ) G H r (Q) x iT" 1 /* (r) : 7 = tr D («)} , 

and note that we have the following dense and compact embeddings V 1 <-» V 2 , for any ri > r2 > 1/2. 
The linear subspace V is densely and compactly embedded into Y, for any r > 1/2. We emphasize that 
V is not a product space and that, due to the boundedness of the trace operator trD, the space V r is 
topologically isomorphic to H r (f2) in the obvious way. Thus, we can identify each u G H r (f2) with a 
pair ( u, trD (u)) G V. Finally, note that both spaces H r (Q) and V are normed spaces with equivalent 
norms. 

The following definition of weak solution to problem (|1.8I) - (|1.10[) is taken from [53] (see, e.g., [2"H 
Definition 2.1], for the more general case). 

Definition 2.1. Let T > and (w ,7o) G Y = L 2 (Vi) x i 2 (r). TTie pair Q{t) = (it(t), 7 (t)) is said to &e 
a (global) weak solution 0/ (03)-(0I23^ 071 [0> ^1 '/> / or almost all t G (0, T], 7(1) = itir(i)) a^d C fulfills 

( G C([0,T];y)nL 2 (0,T;V 1 ), 

a t C G L 2 (0,T;(V 1 )*), ueHl oc ((0,T];L 2 (Q)), 

7 g fri c ((o,r];i 2 (r)), 

such that the following identity holds, for almost all t G [0, T] and for all £ = (x, VO G V , 

(ftC,Ocvi)-,vi + v x > + (f(u), x ) + M)mv) = o. (2.4) 

Moreover, 

C(0) = (m ,7o) =: Co a-e. m Y. 
The map C = (w, 7 ) is a weafc solution on [0, oo) (i.e. a global weak solution,) if it is a weak solution on 
[0,T], for allT> 0. 

Remark 2.2. It is important to observe, for the weak solutions of Definition 12. 1\ that 7 o = uir (0) need 
not be the trace of uq = u\q (0) at the boundary, and so in this context the boundary equation (|1.9[) 
is interpreted as an additional parabolic equation, now acting on the boundary T. However, the weak 
solution does fulfill j(t) = trr>u(t), for almost all t > 0. 
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The existence part of the (global) weak solutions is from [29] Theorem 2.6], and the continuous de- 
pendence with respect to the initial data Co, local Lipschitz continuity on Y, uniformly in t on compact 
intervals, and the uniqueness of the weak solutions follow from [24J Proposition 2.8] (cf. also [29J Lemma 
2.7]). 

Theorem 2.3. Assume l{1.4\ ), \ 1. 5\) and {1. 7| ) hold. For each Co = (uo,7o) G Y, there exists a unique 
global weak solution in the sense of Definition \2.1\ Moreover, the following estimate holds, for all t > 0, 

HCWIly + f t +l \\C{s)\\l^s < CHColl^e-^ + C, (2.5) 

for some positive constants p, C > 0. Furthermore, let C(i) — (u(t),j(t)) and 8(t) = (x(f), ip{t)) denote 
the corresponding weak solutions with initial data Co = («o,7o) an d $o = (xoj^o); respectively. Then, for 
all t > 0, 

\\C(t) - d(t)\\ Y <Ce et \\Co-0 o \\ Y , (2.6) 
where C = C (i?) > is such that \\( \\ Y < R, \\6 \\ Y < It- 
Proof. Since the proofs in [25], [21], [25] involve quite different assumptions on the nonlinearity other 
than the ones in the statement of the theorem, we will sketch a short proof of (|2.5p . This is the main 
estimate on which the proof for the existence of a weak solution is based on (of course, (|2.5p can be 
rigorously justified using a suitable Galcrkin discretization scheme). To this end, testing (|2.4[) with C, 
and appealing to (|2.1I) . we deduce the following inequality 

™ IK Wily + (l - ^) (||V« (t)\\ 2 + \\u (t)\\ 2 L2(r) ) < C, (2.7) 

for all t > 0, where we recall that /i G (0, A] . Exploiting now the continuous embedding V 1 °-> Y, (|2.5I) 
follows from the application of Gronwall's inequality (see Proposition [5T] Appendix) to (|2.7p . The claim 
is proven. □ 

Remark 2.4. Theorem (|2.3p still holds if we keep (|1.7|) and we drop the assumptions Q1.4p . (|1.5p . and 
replace them by the following: 

m\y\ P -C f <f(y)y< m \y\ p + C f , (2.8) 

for some rji, n 2 > 0, Cf > and any p > 2. In this case, the same weak formulation (|2.4[) must be 
satisfied a.e. on [0,T], for all £ = (y, V) G V 1 , with y G L p (fi) (see, e.g., [35] [H]). Finally, we note 
that without assumption (jl.7|) . the uniqueness of weak solutions (given in Definition 12. 1[) is not known in 
general (see [24]). 

Corollary 2.5. Let £/ie assumptions of Theorem \2.3\ be satisfied. We can define a strongly continuous 
semigroup 

So (t):Y^ Y 

by setting, for all t > 0, 

S (t)Co :=C(*) 

where Q (t) = (u (t) , ttir (0) * s ^ e unique weak solution to problem il.8\) - ST.10\) . 

The existence of a bounded absorbing set in V 1 was shown for the first time in j25J Theorem 2.8] and 
the existence of the global attractor for (|1.8|) - (|1.10|) can be found in [24], [25]. The following theorem 
concerns the existence and regularity of the global attractor Aq admitted by the semiflow So and is taken 
from (25] Theorem 2.3]. The proof relies on a uniform estimate which states that problem (|1.8[) - (|1.10[) 
possesses the Y — V 2 smoothing property and exploits (|2.5[) . 

Theorem 2.6. The semiflow So possesses a connected global attractor Ao in Y , which is a bounded subset 
of V 2 . The global attractor Ao contains only strong solutions. Finally, So also admits an exponential 
attractor A^o which is bounded V 2 and compact in Y. 

Remark 2.7. The boundedness of Ao in V 2 , shown in [25] Theorem 2.3], is essential for the proof of 
the continuity property at e = of the global attractors associated with problem (|1.1[) - (|1.3[) . The last 
assertion follows from results in |28[ Theorem 4.2], where (|1.8[) - (|1.10[) is a special case of a phase-field 
system endowed with dynamic boundary conditions. 
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3. The hyperbolic relaxation problem 

In this section, we study the hyperbolic relaxation problem (jl,l[) - (|1.3[l with e £ (0, 1]. Our first goal 
is to prove the existence of a global attractor for (|l.l| - (|1.3j) . As indicated in [51], semigroup methods 
are applied to obtain local mild solutions whereby a suitable estimate is used to extend the solution to 
a global one. We will offer a detailed presentation on the well-posedness of the hyperbolic relaxation 
problem in this section for the reader's convenience. The solution operators define a semiflow on the 
phase space and because of the continuous dependence estimate on the solutions, the semiflow is locally 
Lipschitz continuous, uniformly in t on compact intervals. Further estimates are used to establish the 
existence of an absorbing set for the semiflow. As discussed above, we will follow the decomposition 
method in [33] to obtain the existence of the global attractor in _ff 1 (f2) x L 2 (tt) for the corresponding 
semiflow S E , for each e £ (0, 1]. The (optimal) regularity result for the global attractors A e and a proof 
of their continuity properties conclude the section. 

3.1. The functional framework. Here, we consider the functional setup associated with problem (|1.1|> - 
(|1.3p . The finite energy phase space for the hyperbolic relaxation problem is the space 

He = ff^O) x L 2 {Vt). 

The space % £ is Hilbcrt when endowed with the e-wcighted norm whose square is given by, for ip = 
(u,v) £H £ = H 1 ^) x L 2 (n), 




As introduced in [51] (cf. also [16]), Ar : L 2 {Vl) — > L 2 (fl) is the Robin-Laplacian operator with domain 

D(A R ) = {u 6 H 2 (fl) : d n u + u = on T}. 

Easy calculations show that the operator Ar is self-adjoint and positive. The Robin-Laplacian is extended 
to a continuous operator Ar : H 1 ^) —> (iJ 1 (il)) , defined by, for all v £ H 1 ^), 

(-A R U,W) = (Vlt, Vw) + (U,V) L 2(T). 

Next, [HI EH also define the Robin map R : H S (T) -> H s+ ^(n) by 

Rp = q if and only if Aq = in f2, and d n q + q = p on F. 
The adjoint of the Robin map satisfies, for all v £ 1 (O) , 

R* Arv = —v on r. 
Define the closed subspace of H 2 (tt) x i7 1 (il), 

V e := {(u, v) £ H 2 (fl) x H 1 ^) : d n u + u = -v on T}. 
endowed with norm whose square is given by, for all (p = (u, v) £ T> e , 

IMII> e : = Il u ll2 + Iklli- 

Let D(A e ) = T> e (note that £-dependance does not enter through the norm of 2? e , but rather in the 
definition of A e below). Define the linear unbounded operator A e : D{A S ) — > % £ by 

Ae '- = {lA R I(A R i?tr D -l)) ' 

where trD denotes the Dirichlet trace operator (i.e., trD(v) = w|r)- Notice that if (u,v) £ T> e , then 
u + RtTT>(v) £ D(Ar). By the Lumer-Phillips theorem (cf., e.g., [45j Theorem 1.4.3]) and the Lax- 
Milgram theorem, it is not hard to see that, for all e £ (0, 1], the operator A e , with domain T> e , is an 
infinitesimal generator of a strongly continuous semigroup of contractions on *H E , denoted e Aet . 
Define the map T : H £ —> H £ by 

for all cp = (u,v) £ H £ . Since / : 1 (SI) — > L 2 (il) is locally Lipschitz continuous [SU cf., e.g., Theorem 
2.7.13], it follows that the map IF : H e — >• ~R e is as well. 
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The hyperbolic relaxation problem (jl.l|) - (|1.3[) may be put into the abstract form in T-l e , for tp(t) 
(u(t),u t (t))^, 



(1 



Lemma 3.1. For each s £ (0, 1], the adjoint of A z , denoted A*, is given by 

A* ' 1 

with domain 



±A R -i(A R i? tr D - lW ' 



D(A* S ) := {( X) V) G # 2 (^) x H^Sl) :d nX + X = -r/> on T}. 
Proof. The proof is a calculation similar to, e.g., [H Lemma 3.1]. □ 

3.2. Well-posedness for the hyperbolic relaxation problem. The notion of weak solution to prob- 
lem pJ|) - (TO)) is as follows (see, 

Definition 3.2. A function ip = (u,ut) : [0,T] — > T~i € is a weak solution of \3.1\) on [0, T], if and only if 
J r (y(-)) £ L 1 (0,T;'H e ) and ip satisfies the variation of constants formula, for all t £ [0,T], 

<p(t) = e A ^ Vo + f e A ^~ s *>T{ip{s))ds. 



It can be easily shown that the notion of weak solution given in Definition ^. 21 is also equivalent to the 
following notion of a weak solution (see, e.g., [SJ Definition 3.1 and Proposition 3.5]). 

Definition 3.3. Let T > and (uo,ui) £ T~L e - A map ip ~ (u,ut) £ C([0,T];W e ) is a weak solution of 
H3.1\) on [0,T], if for each 9 = (XjVO G D{A* e ) the map t h-> (ip(t),9)-u e is absolutely continuous on [0,T] 
and satisfies, for almost all t £ [0,T], 

±(p(t),6) ne = (<p(t),A* e 9) He + {T{fP{t)),&)u.. (3.2) 

The map ip = (u,Ut) is a weak solution on [0,oo) (i.e., a global weak solution,) if it is a weak solution on 
[0,T], for allT> 0. 

The above definitions are equivalent to the to the standard concept of a weak (distributional) solution 
to (fl~il- (jl~3l) . 

Definition 3.4. Let e 6 (0,1]. A function ip = (u,u t ) : [0,T] — > "H e is a weak solution of h3.1}) (and, 
thus of fPU-iTO)) on [0,T], if 

ip = («,ut) £ C([0,T];H £ ), u t £ L 2 ([0,T] x V), 

and, for each ip £ H 1 (n) , (u t , ip) £ C 1 ([0, T]) with 

- (eu t (t) , V) + (Vu (t) , V^) + (u t (t) , V) + (ut (0 + u (t) , $ L2(r) =-</(« (*)) , ^) , (3.3) 
/or almost all t G [0, T] . 

Indeed, by [H Lemma 3.3] we have that / : H 1 (fl) — > L 2 (0) is sequentially weakly continuous and 
continuous, on account of the assumptions (|1.4p - (|1.5[) . Moreover, (ip t , 9) £ C 1 ([0,T]) for all 9 £ D(A*), 
and (|3.2j) is satisfied. The assertion in Definition 13.41 follows then from the explicit characterization of 
D{A*) and from Proposition 3.4]. 

Finally, the notion of strong solution to problem (|l.l|) - ([1.3j) is as follows. 

Definition 3.5. Let (po = (uo,ux) £ T> s , e > 0, i.e., (uq,Ui) £ H 2 (n) x H 1 ^) such that it satisfies the 
compatibility condition 

d n uo + uq + u\ = 0, on r. 
A function ip (t) = (u (t) , u t (t)) is called a (global) strong solution if it is a weak solution in the sense of 
Definition \3.4[ and if it satisfies the following regularity properties: 

ip £ L°°(0,c5o;2? £ ), <p t £ L°°(0,oo;?4), 
U tt G L°°(0,oo;L 2 (f>)), u tt € L 2 (0, oo; L 2 (P)) 



(3.4) 
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Therefore, (p (t) = (u(t),u t (t)) satisfies the equations ftl.l\) - [T73\) almost everywhere, i.e., is a strong 
solution. 

We can now state the main theorems of this section. 

Theorem 3.6. Assume \1.4\ ) and il.5\) hold. For each e G (0, 1] and ipo = (uq,ui) G H £ , there exists a 
unique global weak solution (p = (u,u t ) G C([0, oo); H. £ ) to U. 1\) - H77S\) . In addition, 

d n ue Lf oc ([0,oo) x T) andu t G if oc ([0,oo) x T). (3.5) 

For each weak solution, the map 

*^M*)llL+ 2 / F(u(t))dx (3.6) 
Jn 

is C 1 ([0,oo)) and the energy equation 

^Ah(t)\\n e + 2 f Q F(u(t))dx} =-2\\u t (t)\\ 2 -2\\u t (t)\\ 2 L2{r) (3.7) 

holds (in the sense of distributions) a.e. on [0,oo). Furthermore, let ip(t) = (u(t),u t (t)) and 9(t) = 
(v(t),Vt(t)) denote the corresponding weak solution with initial data (po = (uq,ui) G ~H £ an d $o = 
(wo,«i) G We, respectively, such that \\tpo\\u < R, \\9o\\-h — F. Then there exists a constant v\ = 
v\{R) > 0, such that, for all t > 0, 

\Ht) - 9(t)f ne + £ (Hi** (r) - vt (r) || 2 + \\u t (r) v t (r) ||| 2(r) ) dr (3.8) 

< e^yo - e a \\ 2 n! . 

Theorem 3.7. For each e G (0, 1] and (uq,ui) G T> £ , problem il. 1\) - H77S\) possesses a unique global strong 
solution in the sense of Definition \3.5l 

Remark 3.8. The proof of Theorem 13.71 is outlined in [51] (cf., also [16]) when e = 1. 

Proof of Theorem \3.6[ We only give a sketch of the proof. 

Step 1. As discussed in the previous section, for each e G (0,1], the operator A £ with domain D(A £ ) = 
T> £ is an infinitesimal generator of a strongly continuous semigroup of contractions on Ti. £ , and the map 
T : H £ — > H £ is locally Lipschitz continuous. Therefore, by [52] Theorem 2.5.4], for any e G (0, 1] and for 
any ip = (uo,u\) G H e , there is a T* = T* (\\ipo\\n e ) > 0, such that the abstract problem (|3.1[) admits a 
unique local weak solution on [0, T*) satisfying 

ip e C([0,T*);H £ ). 

The next step is to show that T*(||<y9o||"H e ) = °°- Since the map (|3.6I) is absolutely continuous on 
[0, T*) (cf., e.g., [5] Theorem 3.1]), then integration of the energy equation (|3.7p over (0,t) yields, for all 

te[o,T*), 

r /•* /■* 

2 



11^)11^+2/^(^ + 2/ \\u f (r)\\ 2 dr + 2 ||ut(r)||i a(r) dT 

Jn Jo Jo ^3 

= WMh. + 2 I F(uo)dx. 
Jn 

Applying inequality (|2.2|) to (|3.9j) and applying (|2.3|) to the integral on the right hand side, we find that 
there is a function Q(||^o||w e ) > 0, such that, for all t G [0, J 1 *), 

\Ht)\\ Ue < Q(|MkJ- (3.10) 

Since the bound on the right hand side of p,10[) is independent of t G [0, T*), T*(||(^o||hJ can be extended 
indefinitely, and therefore, for each e G (0, 1], we have that T*(||(^ ||« e ) = oo. 

We now show the boundary property (|3.5p . Applying (|2.2p , (|2.3|) and (|3.10p to identity (|3.9p , we obtain 
a bound of the form, for all (po G H £ and t > 0, in which 

ft 

ht(r)\\l H r)<ir<Q(\\ip \\n e )- 

It follows that u t G L 2 oc ([0,oo) x r). By the trace theorem, u G L°°(0, oo; i7 J (0)) ^ L°°(0, oo; L 2 (r)), 
so u G L 2 oc ([Q, oo) x T). Comparison in (jl.2p yields that d n u G Lf oc ([0, oo) x T). 
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Step 2. To show that the continuous dependence estimate (|3.8[) holds, consider the difference z(t) := 
u(t) — v(t), t > 0. We easily get 

±\\(z, Zt )||^ + 2 ||z t || 2 + 2||*||| 3(r) = 2{f(v) /(«), z t ). (3.11) 

Since / : if 1 (0) — > L 2 (f2) is locally Lipschitz continuous, then 

2\(f(v)-f(u),z t )\<Q(R)\\z\\ 2 + \\z t \\ 2 , (3.12) 

where R > is such that ||<£>o||-h < H^oll-^ < Combining (|3.11l) and p,12[) produces, for almost 
all i > 0, 

^ t \\{z,z t )f Hs < Q(R)\\(z,z t )\\ 2 He . (3.13) 

Hence, p. 81) follows immediately from (|3.13l) . using the standard Gronwall lemma. This completes the 
proof of the theorem. □ 

In view of Theorem 13.61 the following is immediate. 

Corollary 3.9. Let the assumptions of Theorem \3.6\ be satisfied. Then, for each each e S (0, 1] we can 
define a strongly continuous semigroup 

S £ (t) :H £ ^H £ , 

by setting, for all t > 0, 

S £ (t) (po = cp(t) = (u (t),u t (*)), 
where (p (t) is the unique weak solution to problem (1. HT73\ ). 

3.3. The global attractor A e in H £ . In this section, we aim to show the existence of a global attractor, 
and prove some additional regularity properties. We point out that all the computations we will perform 
below can be rigorously justified by means of an approximation procedure which relies upon the result 
in Theorem 13.71 Indeed, one shall use the usual procedure of approximating weak solutions by strong 
solutions, and then pass to the limit by using density theorems in the final estimates (see, also, |16j). 
Thus, in what follows we can proceed formally. 

We begin our analysis with a uniform estimate for the weak solutions of Theorem 13.61 The estimate 
provides the existence of a bounded absorbing set B £ C W e , for the semiflow S e , for each e € (0, 1]. 

Lemma 3.10. For all e G (0,1] and ipo = (uq,ui) G H £ , there exist a positive function Q, constants 
LdQ > 0, Pq > 0, all independent of e, such that ip(t) satisfies, for all t>0, 

\\f(t)\\ 2 n e <Q(\\<Po\\n e )e-^ + P . (3.14) 

Consequently, the ball B e in 7-L E , 

B e :={<peH £ :\\<p\\n B <Po + l} (3.15) 
is a bounded absorbing set in 'H £ for the dynamical system (S £ (t) ,"H £ ) . 

Proof. Let e £ (0,1] and tpo = (uq,u\) <G H s . For a > yet to be chosen, multiply (|1.1[) by au in 
L 2 (Q). Adding the result to the energy equation (|3.7p above yields the differential identity, which holds 
for almost all i > 0, 

^ jlM&. +ae{ut,u) +2 j^F(u) dzj + 

+ (2 - ea)\\u t \\ 2 + a(u t ,u) + a||u|| 2 + (3 ' 16) 
+ 2 IMIz, 2 (r) +ol{u u u) L 2 [T) +a(f(u),u) =0. 
For each e E (0, 1], define the functional, E e : H £ — > R, by 

E £ (<p (*)) = y (t) || 2 ie + as(u t (t) , u (t)) +2 f F(u (t))dx. (3.17) 

Jn 

It is not hard to see that the map 1 E £ (ip(t)) is C 1 ([0,oo)); this essentially follows from equation (|3 . 6f) 
of Theorem 13.61 First, we estimate, for all r\ > 0, 

a\(u t ,u) L 2 (r) \ < a?]\\u t \\l 2(r) + ^\\u\\ 2 L 2 (r) , (3.18) 
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and with (|2.1|) , we have, 

a|(/W ! u)|>-^-^NI?-aC. (3.19) 
Combining (|3~T^|) with ([335]) - ([339]) gives 

^-E E + (2-a)e\\ut\\ 2 + a(u t ,u) (3.20) 
di 



1 A-/x 

all-- — 

4rj A 

< aC 



j |HI? + (2-cw 7 )||u t ||i a(r) 



Hence, for any r\ > ^ an d an y < a < min{2, -}, then 2 — 77 > and 2 — ar\ > 0, 

w ° :=min { 2 - a ' a (^-^)} >0 ' 
and estimate p.20[) becomes, for almost all t > 0, 

^S E +cJo^ e + (2-a?7)||u t ||| 2(r) < C a . (3.21) 

Applying Gronwall's inequality (see, e.g., [HI Lemma 5]; cf. also Proposition 15. 1[ Appendix) to (|3.21l) 
produces, for all t > 0, 

E e (<p(t)) <^(v(0))e-^ t + C. (3.22) 
We now apply (|2.2[) to (|3.17p to attain the bound, 

S E M > e (i-2) W » + (i_|L_ ^) H? - C. (3.23) 
After updating the smallness condition on a to < a < min{2, ^, 2/i}, we see that for 

• f 1 a 1 Q A ~ M I n 

a* := mln |1 - _, 1 - _ > 0, 

then, for all i > 0, 

^( V (t))> Wl || V (t)||^-G (3.24) 
On the other hand, by estimating in a similar fashion, using (|2.2[) , there holds for all t > 0, 

^(^))<Q(||v(t)ll«.). (3-25) 

Thus, estimate Q3.14p follows now from (|3.24p . (|3.25p and (|3.22p . The assertion (|3.15p is an immediate 
consequence of (|3.14p . This concludes the proof. □ 

Remark 3.11. The following bounds are an immediate consequence of estimate (|3.14p : 

limsup ||¥>(*)|||i s < Po (3.26) 

t— yoo 

and 

(||«t (r)|| 2 + \\ut (r)|£ a(r) ) dr < Q (\\M\nJ ■ (3-27) 

The last bound is found by integrating the energy equation Q3.7j) with respect to t over (0, oo) and 
estimating the result with ([L~gj) . ([T5]) . (|2T2"j) . ([2~3]) and (|3~2"6"|) . 

Remark 3.12. Note that the last assumption ()1.7[) (which is, / (s) > —0, for all s £ M) is nowhere needed 
in the proofs of Theorem 13. 6[ Theorem 13.71 (cf. [51] Theorem 1.1 and Lemma 2.2]) and Lemma [3.101 
It will only become important later (see (|3.28p ) when we establish the optimal regularity of the global 
attractor for the hyperbolic problem (ll.ll) - (|1.3p . 
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The semiflow S e admits a bounded absorbing set B £ in H £ . To obtain a global attractor, it suffices 
to prove that the semiflow admits a decomposition into the sum of two operators, S £ = Z £ + K £l where 
Z £ = (Z £ (t)) t >o and K £ = (K £ (t)) t >o are not necessarily semiflows, but, operators that are uniformly 
decaying to zero, and uniformly compact for large t, respectively. To obtain the compactness property 
for the operator K e , recall that, when fractional powers of the Laplacian arc well-defined, one usually 
multiplies the PDE by the solution and a suitable fractional power of the Laplacian; i.e., (— A) s u for 
some s > 0, then estimates using a stronger norm while keeping in mind the uniform bound on u and 
the null initial conditions. However, in our case, the dynamic boundary condition does not allow us to 
proceed with the usual argument to obtain the relative compactness of K £ . This is because the Laplacian 
equipped with the dynamic boundary condition (|1.2p is not self-adjoint nor positive. In turn, we cannot 
apply the standard spectral theory to define fractional powers of the Laplacian. So to obtain the relative 
compactness of K e , we follow the approach in |44j . The main tool is to differentiate the equations with 
respect to time, and obtain uniform estimates for the new equations. Such strategies also proved useful 
when dealing with a damped wave equation with acoustic boundary conditions [23], or a wave equation 
with a nonlinear dynamic boundary condition |16| , and hyperbolic relaxation of a Cahn-Hilliard equation 
with dynamic boundary conditions [14] . [27j . 

Following an approach similar to the one taken in the above references, first define 

iP(s):=f(s)+ps (3.28) 

for some constant ft > d to be determined later (in this case, ip'(s) > thanks to assumption (JT77J)). Set 
^(s) := J Tp(a)da. Let ipo = (uq.ui) G TL £ . Then rewrite the hyperbolic relaxation problem into the 
system of equations in v and w, where v + w = u. 

ev tt + v t — Av + ip(u) — i'(w) = in (0, oo) x f2, 

d n v + v + v t = on (0, oo) x T, (3.29) 

v(0) = uo, «t(0) = ui + / (0) - f3u in Q, 

and 

ewtt + Wt — Aw + ip(w) = /3u in (0, oo) x fl, 

d n w + w + w t = on (0, oo) x T, (3.30) 

w(0) = 0, u*(0) = -/ (0) + /3m in n. 
In view of Lemmas 13. 131 and 13. 151 below, we define the one-parameter family of maps, K £ {t) : H s — > H e> 

by 

K e {t)<p := (w(t),w t (t)) , 

where (w,wt) is a solution of (|3.30[) . With such w, we may define a second function (v,v t ) as the 
solution of p.29p . Through the dependence of v on w and <^o = (wo,wi), the solution of (|3.29[) defines a 
one-parameter family of maps, Z £ (t) : H. £ — > T-L £l defined by 

Z e {t)<Po ■= («(*),««(*)). 

Notice that if v and w are solutions to p.29|) and p.30p . respectively, then the function u := v + w is a 
solution to the original hyperbolic relaxation problem (|1.1[) - (|1.3[) . 

The first lemma shows that the operators K £ are bounded in 'H £ , uniformly with respect to e. The 
result essentially follows from the existence of a bounded absorbing set B £ in H. £ for S £ (recall (|3.26[1 ). 

Lemma 3.13. Assume 11. 5\) and jl.7\ ) hold. For each e 6 (0,1] and ipo = (uo,Ui) £ H £ , there 

exists a unique global weak solution (w,w t ) £ C([0, oo); H £ ) to problem \3. 30\) satisfying 

d n w G Lf oc ([0,oo) x F) and w t G Lf oc ([0, oo) x F). (3.31) 

Moreover, for all ipo G T~L £ with ||(/?o||% < R for all e G (0, 1], there holds for all t > 0, 

\\K s (t)cp \\n e < Q(R). (3.32) 

The following result will be useful later on. 

Lemma 3.14. For all e G (0, 1] and r\ > 0, there is a function Q v (•) ~ rj _1 , such that for every < s < t 
and tpo = (uo,Ui) G B £ , 

f {\\w t {r)\\ 2 + \\u t (r)\\ 2 ) dr < ^(t s) + Q n (R) , (3.33) 
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where R > is such that Hvoll-w < R, for all £ G (0, 1]. 

Proof. Let e € (0, 1] and ipo = (ito,wi) € % £ , with ||^o||^ < -R- Adding the identity 



-20^{u,w) = -2p{u t ,w)-2P{u,w t ) 



to equation 

d_ 

dt 

produces, for almost all t > 0, 
d 



\\(w,w t )\\i lt +2 / *( W )da;^+2|K||' j + |K||i a(r) =2 J 8(«,«; t > (3.34) 



^IK^)llL+2 ^*Hdo:-2/3( W , U ;)|+2||«; t || 2 + || Wt ||i 2{r) 

= -20(ut,w)- 

Using (|3.32|l , we estimate, for all r/ > 0, 

2y9|<t* t ,to>| <r/ + Q^(fl)||u 4 || 2 . (3.36) 
For each e G (0, 1], define the functional W £ : if^fi) x H^fl) x i 2 (fi) M, 

W £ (i) :=!|( W (t), W4 (t))|| 2 +2 / *(w(t))da ! -2 J 8(ti(t) > i£;(t)). 



Because of ([23]) , ([Tip . (JT3J), (j3~^| . ([3361) and ([332]) , we can easily check that for all £ > and e € (0, 1], 

\W £ (t)\ < Q (R) . (3.37) 
We now combine (|3.35p and p.36[) together as, for all rj > and for almost alH > 0, 

^-W £ + 2\\w t \\ 2 + K||| 2(r) + 2||u t || 2 < r) + (Q n (R) + 2) \\u t \\ 2 . (3.38) 

Integrating (|3.38p over (0, t), and recalling (|3.37p . (|3.27p . gives the desired estimate in (|3.33p . This proves 
the claim. □ 

The next result shows that the operators Z £ are uniformly decaying to zero in Ti. £ . 

Lemma 3.15. Assume il-5) ) and Ji.Tj ) hold. For each e G (0,1] and ipo = (uq,ui) G H £ , there 

exists a unique global weak solution (v,vt) G C([0,oo);H e ) to problem H3.29\) satisfying 

d n v G L 2 oc ([0,oo) x T) and v t G Li 2 oc ([0,oo) x T). (3.39) 

Moreover, for all ipo G T> £ with \\tpo\\ft < R f or a ^ e G (0, 1], there exists uj > 0, independent of e, such 
that, for all t > 0, 

\\Z £ {t)Mu e < Q(R)e~ ut . (3.40) 



Proof. In a similar fashion to the arguments in Section 13.21 the existence of a global weak solution as 
well as (|3.39p can be found. Because of (|3.26p and (|3.32p . we know that the functions (u(t),Ut(t)) and 
(w(t),Wt(t)) are uniformly bounded in H £ with respect to t and e. It remains to show that (|3.40p holds. 
Let e G (0, 1] and cpo = (uq, ui) G H £ , with R > such that ||<^o||^ < R- Observe that, 

2(V(u) - il>(w), v t ) = 1 {2{i){u) - i>{w),v) - ltf(v)v, v)}- 

- 2{(ip'(u) - ip'(w))wt,v) + (ip"(u)u t ,v 2 ). 

Multiply the first equation of p. 291) by 2vt + ctv in L 2 (il), for a > to be chosen later. We find that, 
with ([gilt , 

^ {e||f t || 2 + ae(vt,v) + \\v\\ 2 + 2(ip(u) - i/)(w),v) - (i>'(u)v, v)} + 

+ (2 -ae)\\v t \\ 2 +a(v t ,v) + a\\v\\ 2 + 2||w t || 2 2(r) + a(v t ,v) L 2 {r) + ( 3 - 42 ) 
+ a{ip(u) - ip(w), v) = 2((ip'(u) ~ tp'(w))w t ,v) - (ip"(u)u t , v 2 ). 
For each e G (0, 1], define the functional 

V e : H x (Sl) x H 1 ^) x H^fl) x L 2 (0) -> M, 
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IS 



by 

V&) ■.=e\\v t (t)\\ 2 + ae(v t (t),v(t)) + \\v(t)\\ 2 + 

+ 2(i>( U (t)) ~ ^(w(t)),v(t)) - (ip'(u(t))v(t),v(t)). 

As with the functional E E above, the map t H> V e (t) is ^4C(R>o; K>o). We now will show that, given 
(u,u t ),(w,w t ) € T-L E are uniformly bounded with respect to t and e, there are constants, G\,C<x > 0, 
independent of t and e (possibly depending on R > 0), in which, for all (v, v t ) £ W E , 

Ci||(«,«t)||«.<V e <C 2 ||(t;,« t )||^. (3.43) 

We begin by estimating the products in V E that involve tp; with (11.4p . (|1.5[) . the embedding iJ 1 (SI) 
i 6 (0) and (f3T26|) . there holds 

|(V/(^>l<Cn(i + !MI?)IMIilMI 

1 (3-44) 
< 2ll«ll? + W)ll«H a - 

From assumption (| 1 . T[) and the definition of t/j, cf. Q3.28p . 

2{i>{u) - i>(w), v) > 2(/3 - i?)|M| 2 . (3.45) 

Hence, for f3 sufficiently large, /3 > (C (R) + 2i9) /2, the combination of (|3.44[) and (|3.45j) produces, 

- *(»),«) - (•>'(<•)»,») > 208- - - C(B)||«|| 2 

Then we attain the lower bound on V E , 



So for < a < min{2, A}, set 



-f).M- + (i-£ 



a 1 a . 

ujo '■= mm < 1 , > > 0. 

2' 2 2A ' 



then, for alH > 0, we have that 

V e (t)>u> 2 \\(v(t),v t (t))\\n B - ( 3 - 46 ) 
Now by the (local) Lipschitz continuity of /, and the uniform bounds on u and w, it is easy to check that 

2(V(u) - 1>(w),v) < 2\\1>(u) - 1>(w)\\\\v\\ < Q(R)\\v\\l 

Also, using (jl.4[) . (| 1 . 5|) and the bound (|3.26p . there also holds 

\W(u)v,v)\<Q(R)\\vf v (3.47) 

Thus, the assertion in (|3.43[) holds. Exploiting the fact that 

a\(v t ,v) L 2 (r) \ < ^\\v t \\ 2 L2{r) + ^\\v\\ 2 L 2 (r) , 

we see that (|3.42[) becomes 

— V E + (2 - a)e\\v t \\ 2 + a{v u v) + a\\Vv\\ 2 + |||«||i a(r ) + 



+ (2 - ^)\\vt\\h {r) + a(^(u) ~ VH, v) - W(u)v,v) ( 3 - 48 ) 
< —{tp'(u)v, v) + 2((V>'(u) - ip'(w))w t ,v) - (4>"(u)u t ,v 2 ). 
Recall that < a < min{2, A}, so when we set 

co 3 := min |2 - a, 1, ^ j > 0, 

we write (|3.48p as 

^V E + u a V e < -W{u)v,v)+2((i>'{u)-i>'(w))w t ,v) - {iP"{u)u u v 2 ). (3.49) 
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Using the uniform bound on u and w (recall assumptions (jl.4p . (|1.5p . (|3 . 26[) and (|3.32[) ). there is a positive 
function Q V (R) > 0, depending on 77, such that, for all n > 0, 

- ^' H) u*,«)| < Co (1 + IklK + Hli) H|| Nli (3-50) 

<|NI? + Q,(fl)IKI|V £ . 

The last inequality in the above estimate follows from (|3.46[) . In a similar fashion we estimate using 
assumption (|1.4[) and the bound p.26|) . 

\^"(u)u u v 2 )\ <Cn(l + ||w|li)ll«tlllHli (3-51) 

<f NI? + Q,(£)lkll 2 K. 

Applying ([3^44)1 to (|3~49|) and inserting (|330)) and ([3~5Tj) into (|3~49)l . we then have 

+ w 3 V e - »?||«||i < Qr, (R) {\\u t \\ 2 + \\w t \\ 2 ) V e . (3.52) 
There is a sufficiently small 77, precisely, < r\ < OJ3/2, so that (|3.52[) becomes 

^V £ + rjV s < Q n (R) (\\u t f + \\w t f) V £ . (3.53) 

At this point, we remind the reader of Lemma 13.141 Applying a suitable Gronwall type inequality (see, 
e.g., jUJ Lemma 5]; cf. also Proposition 15. 1[ Appendix) to (|3.53|) yields 

V s (t) < K(0)^ (Ji) e^ /2 . (3.54) 

By virtue of (j3~43)) . for all s G (0, 1], 

V e (0) < Q(R)\\(v(0),v t (0))\\l c 

< Q(R)(\\u Q \\ 2 1 +e\\u 1 +f(0)-(5u \\ 2 ) 

< Q (R) , 

for some positive function Q independent of e. Therefore (|3.54[) shows that the operators Z E are uniformly 
decaying to zero. □ 

The following lemma establishes the uniform compactness of the operators K £ . 

Lemma 3.16. For all tp = (u ,ui) G H e such that ||yo||-w < R f or a ^ e G (0? 1]; the following estimate 
holds: 

\\K e {t)<po\\v. < Q(R), 
for all t > 0. Furthermore, the operators K £ are uniformly compact in T-L e . 

Proof. Let e G (0, 1] and let <^o = (uoi u i) G %e with R > such that H^oII-h < -R- Differentiate (|3.30l) 
with respect to t and set h = w t . Then h satisfies the equations 

ehu + h t — Ah + ip'(w)h = f3u t in (0, 00) x Q, 

d n h + h + h t = 0, on (0,oo) xT, (3.55) 

h(0) = w t (0) , ht(0) = w u (0) in fl. 

Note that, by the choice of data in p.30|) . we actually have h (0) = — / (0) + f3u and h t (0) = 0. Multiply 
the first equation of p.55|) by 2h t + ah, where a > is yet to be determined, and integrate over fi. 
Adding the result to the identity 

2(^'(w)h,h t ) = j t W{w)h,h) - (ip"(w)wt,h 2 ) 

produces 

^ {e\\h t \\ 2 + ae(h t , h) + \\h\\l + W(w)h, h)} + 
+ {2-ae)\\h t \\ 2 +a(h tl h) + a\\h\\\+ 
+ 2 ll^lli, 2 (r) + a(ht,h) L z(r) + a(ip'(w)h, h) 
= (ip"{w)w t ,h 2 ) + 2/3(u t ,ht) + a/3(u u wt). 
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For each e E (0, 1], define the functional 

* £ : ff^Q) x H 1 ^) x L 2 {VL) -> R 

by 

tt e (t) := e\\ht(t)\\ 2 + ae(/i t (t), fc(*)> + \\h{t)\\\ + (ij/(w(t))h(t),h(t)). (3.57) 
The map 1 1-> v E , £ (t) is v4C(R>o; K>o)- Because of the bound given in Lemma r3.131 we obtain the estimate 
similar to (|3.47|l . 

a\{ip'(w)h,h)\ < aQ(R)\\h\\j. (3.58) 

Obviously, we have 



as\(h t ,h)\<^\\h t \\ 2 + ^.\\h\\l (3.59) 



After combining (|3.57[) - (|3.59[) . we find 



* £ > ( 1 - 77 ) e\\htV + ( 1 - Y\ - aQ ^ ) 



1 



Hence, when 

< a < min I2, ^ + Q(fl) 
then, 

W4 (i?) := min {l - -, 1 - — - aQ(R)j > 0, 

thus, for all t > 

*e(*) >W4||C»(t),/lt(*))||«.. (3-60) 

On the other hand, again with (|3.58l) , 



* £ < ( 1 + £ ) e||M 2 + f 1 + 7^ + «Q(i?)) ll^ll?, 



2/ " " V 2A 



and with 



w 5 (-R) :=max|l + -,l + — + aQ(R)j, 
an upper-bound for \I> £ is given by, for all t > 0, 

* e (*)<W5||(M*),^(*))ll%- (3-61) 
Using the bounds found in (|3 . 26[) and (|3.32[) , we estimate the following terms from (|3 . 56[) . for all 77 > 0, 

a\(h t ,h) LHr) \ < arj\\h t \\ 2 L2{r) + ^\\h\\ 2 L2{r) , (3.62) 

and 

2f3\(u t7 h t )\+a(3\(u t7 w t )\ < Q(R)\\h t \\ + Q(R) 

„ „9 , , (3-63) 

<il\\h t \\ 2 + Q n {R). V ' 

Also, similar to (|3.51[) . but when we now employ (|3.60[) . we have that, for all 77 > 0, 

W(w)w t ,h 2 ) < Q n (R)\\w t \\* B . (3.64) 

Combine (|3.62l) - (|3.64p with (|3.56p and obtain the following estimate (note that when 2 — a — 77 > 0, we 
have (2 — a — ij)e < 2 — as — rf): 



-tt e + (2 - a- V )s\\h t \\ 2 + a(h t , h) + a\\Vh\\ 2 + a - — J ||/i||£ 2(r) + 

+ (2 - ari)\\h t \\l 3{r) + a(i>'(w)h, h) < Q n (R)\\w t \\* e + Q V (R). 
With some j < 77 < 2 now fixed, then, for 



we have 

elf 



< a < min < 2 - 77, 1, — > and wg :— 1 , 

I V) 477 

d * £ + ^ 6 * £ + (2 - cwjJHMW) ^ G ( E ) + ( 3 - 66 ) 
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An immediate consequence of (|3 . 2T[) is the bound on the following integral 

IK(r)|| 2 dT<Q(i?). 



Applying a suitable version of the Gronwall inequality (see, e.g., [32j Lemma 2.2]; cf. also Proposition 
15.11 Appendix) it follows that 

# e (t) < Q (R) * e (0)e^ 6t/2 + Q(R). (3.67) 

Using ([3~60]) and (j3~6Tjl . and the fact that * e (0) < u) 5 \\(h(0), h t (0))\\^ e < Q (R), we arrive at the bound 

\\w t (t)\\l + 4mt(t)\\ 2 <Q(R), (3.68) 

for all t > 0, e <E (0, 1] and tpo <G T~L €l with R > such that ||yo||^ < R- 

In order to bound || (w, wt)\\v c , we need to bound the term ||io||2- We have owing to standard elliptic 
regularity theory (see, e.g., [33 Theorem II. 5.1]), that 

Kt)||a < C (\\Aw{t)\\ + \\d n w(t)\\H^(r)) ■ (3-69) 

Thus, using the first equation of (pT30)) . the bounds ((3~26|) . (|3~32|) and (|3~68)) . and also (fO]) . ([LB) and 
([3T28]) . we have 

|| Aw (*) || < v/i||itf« (t) || + |K (*) II + \\1>(w (*))ll + P\W (*) II < QW- (3-70) 
Also, by ([3T68]) . we have that w t G L°° (M> , # 1/2 (L)) . Thus, from the second equation of (pOD]) . 

Hft.il> (t) || H i/ 2(r) < ||w (t) || H i/ 2(r) + Ikt W || H i/ a( r) < Q{R)- (3-71) 
Combining ([3~70)) and f3~7T) with ([3T69]) . and also applying ([3~68)l . proves that for all i > 0, 

\\(w(t),w t (t))\\v e < Q(R). 

It follows that the operators K e are uniformly compact (with t c = 0). □ 
Next, we will discuss regularity properties of the weak solutions. 

Theorem 3.17. For each e G (0,1], there exists a closed and bounded subset C £ C T> £ , such that for 
every nonempty bounded subset B C H e , 

dist nc (S e (t)B,C e ) < Q(\\B\\ ne )e-^, (3.72) 

where Q and oj > are independent of e. 

Proof. Let e <G (0, 1]. Define the subset C e of V e by 

C £ :=WeV e : y\\ Vs <Q{R)}, 

where Q(R) > is the function from Lemma 13.161 and R > is such that ||yo||H £ ^ R- Let now 
ipo = (uo, u i) € &s (endowed with the same topology of H e ). Then, for all t > and for all ipo £ B e , 
S £ (t)(po = Z e (t)tpa + K e (t)tpa, where Z e (t) is uniformly and exponentially decaying to zero by Lemma 
13.151 and, by Lemma 13.161 K e (t) is uniformly bounded in T> e . In particular, there holds 

dist ne (S e (t)B e ,C e ) <Q(R)e^ t . 

(Recall, io > is independent of e due to Lemma [3.151) . 

Recall that, by Lemma l3.101 we already know that for each e € (0, 1] and for every nonempty bounded 
subset B of He, 

distn e (Se{t)B,B £ ) <Q{R)e-^\ 

for all t > 0. In light of these estimates, (|3.72[) can now be accomplished by appealing to the transitivity 
property of the exponential attraction (see, e.g., [22j Theorem 5.1]). Note that (|3.72[) entails that C e is a 
compact attracting set in H e for S e (t). The proof is finished. □ 

By standard arguments of the theory of attractors (see, e.g., [3~H [H]), the existence of a compact 
global attractor A e C C e for the semigroup S e (t) follows. 
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Theorem 3.18. For each e £ (0, 1], the semiflow S £ generated by the solutions of the hyperbolic relaxation 
problem 11. 1)) - (T73\) admits a unique global attractor 

A E =u(B E ) := p| \JS E (t)B E * 

s>t t>0 

in Ti. £ . Moreover, the following hold: 

(i) For each t > 0, S £ (t)A £ = A £ , and 

(ii) For every nonempty bounded subset B ofH e , 

lim dist-H £ (S E (t)B,A E ) = 0. (3.73) 

t— >oo 

(Hi) The global attractor A E is bounded in T> £ and trajectories on A E are strong solutions. 

Remark 3.19. We can extend all the results in Sections 13.2113.31 (with the appropriate modifications, see 
[H1I2Z]) t° the case when the linear boundary condition (jl.2[) is replaced by 

d n u + g (u) + ut = 0, on (0, oo) x T, 

such that g £ C 2 (K) satisfies 

Iff" (s) \<C g (l+ \s\ 2 ) , g (s) > -6 g , g(s)s>s 2 - C' g , 
for all sel, and some constants C g > 0, C g > 0. 

3.4. The upper-semicontinuity of A £ for the singularly perturbed problem. This section con- 
tains one of the main results of the paper, the proof of the upper-semicontinuity of the family of global 
attractors given by the model problems for e £ [0, 1]. Recall that in the case e = 0, the limit parabolic 
problem, admits a global attractor Ao that is bounded in V 2 . Naturally, we will study the continuity at 
e = 0. For eG (0, 1], we know that the hyperbolic relaxation problem admits a global attractor A £ in T> £ . 
However, the spaces involved with the parabolic problem invoke the trace of the solution on the bound- 
ary r, whereas the spaces involved with the hyperbolic relaxation problem do not contain prescribed 
traces. Before we lift the global attractor Ao for the parabolic problem into the finite energy phase space 
for the hyperbolic relaxation problem, we need to make an extension of Ti. £ so that it also includes the 
information of the traces of u and Ut- 

To begin, we recall that the natural phase space for the parabolic problem (|1.8[) - (|1.10[) is Y = 
L 2 (il) x L 2 (T), while the finite energy phase space for the hyperbolic relaxation problem (|l.ll) - (|1.3|) 
is H £ = H 1 ^) x L 2 (il). Thus, we need to find a suitable extension of the phase space for the hyperbolic 
relaxation problem so that, when we lift the parabolic problem, both problems will be situated in the 
same framework. A natural way to make this extension is to introduce the space 

x = h x {si) x i 2 (r), 

then the extended phase space for the hyperbolic relaxation problem 

X £ = X x Y = x L 2 (T) x L 2 (fl) x L 2 (T). 

The space X £ is Hilbert when endowed with the e-weighted norm whose square is given by, for all 
C = (u,l,v,S) £ X E , 

l!CllL:-||«ll? + ll7lli 2 (r)+^IH| 2 + s||^ ( r)- 
It is then in the space X £ that we can lift Ao and estimate the HausdorfF scmidistance between (an 
extension of) A £ and CAo (for a proper lifting map £) with the new extended topology. However, it 
must be noted that the lifted attractor CAo is n ot necessarily a global attractor when set in the extended 
phase space. Finally, the topology that we will use to show the convergence of the attractors at e — 
will be defined with the four-component norm of X £ . 

For both problems, we also recall that trajectories on the attractor are strong solutions due to the 
regularity results obtained in Sections [2] and 13.31 (see Theorems 12.61 and I3.17|) . The regularized phase 
space T> £ for the hyperbolic relaxation problem is isomorphically extended to 

V £ := { (u, 7, v, 5) £ H 2 (O) x B 3 ' 2 (r) x H 1 (O) x H 1/2 (V) : 7 = tr D («) , (3.74) 

S = trD (v) , d n u + 7 = —5 on T} . 
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Of course, T> £ C V 2 x V 1 and the injection T> £ X £ is compact. Recall that, for each (uq,u\) £ T> £ , 
problem (|l.ip - (|1.3j) generates a dynamical system (S £ (t) ,T> e ) of strong solutions (cf., Theorem I3.7( see 
also [5T]). By appealing once more to the continuity of the trace map tro : H s (Q) H 8 ^ 1 ! 2 (r), s > 1/2, 
and exploiting the results in Section \3. 21 it is not difficult to realize that we can extend the scmiflow S £ (t) 
to a strongly continuous semigroup 

S £ (t) : V £ -> T> E , (3.75) 

such that S £ (t) is also Lipschitz continuous in T> £ , endowed with the metric topology of V 2 x V 1 (sec 
Lemma [3.201 below). Recall that, by definition for p,q> 1, 

VxV« = t(u,i,v,6) £ H p (fi) x flP-Va (r) x H q (Q) x ff^Va (T) : 7 = tr D («),<* = tr D («)} , 

see Section [5] (as before, V p x V 9 is topologically isomorphic to H p (0) x if 9 (fi)). 

Lemma 3.20. Let <^o,#o G sucft t/iai ||v?o||£S < R> \\@o\\t> ^ ^> / or every e G (0,1]. TVien the 
following estimate holds: 

S e {t)<po-S e (t)0 o _ < QMe^Wipo-OoWz , (3.76) 
where Q (R) > and v\ > are independent of e > 0. 

Proof. Let = (ui(t),uxir (<) ,9t«i(*),8t«x|r(*)) and 6>(i) = (u 2 (i), u 2 | F (i) , d t u 2 (t),d t u 2 \r{t)) denote 
the corresponding strong solutions with initial data (/?o and #0, respectively. Then the difference u (t) := 
u\ (t) — u 2 (t) satisfies 

-Am (t) = f (u 2 (t)) - /' (ui (t)) - u t (t) - eu tt (t) , a.e. in R + x ft, , , 

d n u (t) +u(t) = -u t (t) , a.e. in 1 + x T, [ ' ' 

subject to the initial condition 

u(0)=«i (O)-ua (0). 

Setting v :— dtU\ — dtU2, we have (vt,ip) £ C 1 ([0, T]) for every tp £ H 1 (f2) (see the definition of strong 
solution). Then v solves the following identity 

^ (sv t (t) , yj) + {Vv (t) , W) + (vt (t) , V) + (vt (t) +v{t), i>) L2{r) 

= - (/' («i (*)) ~ /' («2 (*)) ui (t) , 1>) - (/' («a (t)) u (t) , </>) , 
for almost all t £ [0, T] . Testing with tp = Vt, we obtain 

^ {e Ikf + l|V«|| 2 + ||«||" 2(r) } + |h|| 2 + IMI^r) (3.78) 

= - (/' (Ml) - /' («2)«l,ft) - (f (u 2 )u,V t \ ■ 

We can bound the terms on the right-hand side in a standard way, 

(f (til)- / («2)«i,ft) +(/ (u 2 )u,v t \ 

<Q(Noo)Hl a + ^IHI 2 

(which follows easily on the account of the fact that || (iij (i) , dtUi (t))^ < R, i = 1, 2, and the embedding 
if 2 (fi) <-> C° (H)), then insert them into ([3~75|) . By virtue of ([575]) we get 

£|kWH 2 +(||V«(t)|| 2 + ||t;(t)||^ (r) ) (3.79) 

< Q (R) e" 1 ' ll^o - 0o|&. + e (°)l| 2 + (llW) - flolll.) , 
for almost all i G [0, T] . It remains to notice that, from (|3.77[) . there holds for every e £ (0,1], 

e Hvt (0)|| 2 < - e (||Au(0)|| 2 + ||/' (u 2 (0)) - /' («! (0))|| 2 + (0)|| 2 ) (3.80) 

S 11^0 - ^oll-D 
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Summing up, we obtain from (|3.T9[) - (|3.80[) . that 

e \\v t (t)\\ 2 + \\v (t)\\l < Q^e** y Q O ||! . (3.81) 



£ 

We can now bound the term \\ui (t) — u-2 (t) \\2- As before, owing to standard elliptic regularity theory, 
we have in p.77[) . using f|3 .81 [) . that 

\\u{t)\\l < c(£ 2 \\v t (t)f+\\f'(u 2 (t))-f'( Ul (t))f + \\v(t)\\^ (3.82) 

< Q(R) \\<p -6 fv e . 

Finally, (|3~8T|) - (|X82|) together with the fact that the trace map H s (Q) H 8 ' 1 / 2 (T), s > 1/2, is bounded 
yields the desired inequality (|3.76p . □ 

By Lemma [3.201 the family of global attractors {„4 £ } e6( - ^ C T> £ can be naturally extended to the 

family of compact sets i 

Ae = |(u,7,«,5) e V £ : (u,v) E A e \ (3.83) 

which are bounded in T> £ and compact in X e . Note that we do not claim that A e is a global attractor 
for (S e (t) , X e ), see Remark 13 . 2 1 1 below . Also, it is in the space V 2 x Y C X e where we lift the parabolic 
problem. Since the global attractor Ao for (|1.8|) - (|1.10|) is a bounded subset of the space V 2 C C (fi) x C (F) 
(since SI C K 3 ), the canonical extension map 

£ : V 2 -)• Y (3.84) 

is well-defined with 

(u, uir) !->• (Ait — f(u), —d n u — uir)j (3.85) 

and so the corresponding lift map 

£:V 2 4V 2 xF (3.86) 

is defined by 

(u, u\ T ) n- (u, u\ r , Aw - f(u), -d n u - u\ T )- (3.87) 

Let Aq denote the global attractor of the limit parabolic problem (see Fheorem 12. 6|) and let A e , 
e € (0, 1], be the sets defined in (|3.83|) . Define the family of compact sets in X £ by 

A e := { 4° ;= CA ° f ° r £ = ° , (3.88) 
\ A e for £ G (0,1]. V ; 

Remark 3.21. Fhe compact set A e is not a global attractor for S e (t) acting on the phase-space X e since 
traces of functions in L 2 (f2) are not well-defined in L 2 (r). By construction (|3.83|) . A £ is only topologically 
conjugated to the global attractor A e associated with the dynamical system (S £ ,H £ ) ■ 

The main result of this section can be now stated as follows. 

Theorem 3.22. The family {A e } eg j ^ , defined by L8. 88]) . is upper- semicontinuous at £ = in the 
topology of X\ . More precisely, there holds 

lim dist^ (A e , Ao) := lim sup inf \\a — b\\x 1 =0. (3.89) 

Proof. Our proof essentially follows the classical argument in 33,(34] and also [40j Theorem 3.31]. Of 
course, modifications are required to account for the terms on the boundary. Let £ = (u, 7, v, S) G A e 
and ( = (u, 7, v, 6) e Ao. We need to show that 

sup inf (\\u - u\\j + ||7 - 7lli2(r) + 

\ 1/2 

+ ||7;-f;|| 2 + ||<5-5||| 2(r) J ->0 ase^O. 
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Assuming to the contrary that (|3.90[) did not hold, then there exist r] > and sequences (e„)„ e N C (0, 1], 
(Cn)neN = ({u n ,Jn,v n ,S n )) ne n C A Sn , such that e„ -> and for all neN, 

inf _ (K - «||? + ||7» - 7ll| 2( r) + K - "f + \K - S\\h {r) ) > vl (3.91) 

By Theorem 13.171 the compact sets A Erl are bounded in the space T>\ (see (|3.74[) with e = 1) and we have 
the following uniform bound, for some positive constant C > independent of n, 

ll M "ll2 + Il7«ll#3/2( r ) + H^nlll + ||£re||fl-l/2(r) < C. 

This means that there is a weakly converging subsequence of (Cn)neN (not relabelled) that converges to 
some (u* , 7*, v* , 6*) weakly in T>\, By the compactness of the embedding T>\ ^ X\, the subsequence 
converges strongly in X\. It now suffices to show that (u*, 7*, v*, 5*) £ Ao, since this is a contradiction 
to (EH]). 

With each („ — [u n , 7«, v n , 6 n ) £ A 6n , then, for each n £ N, there is a complete orbit 
(u n (i),uj^(i),u?(i),u^ r (i)) teR = (Se n (t)(u„,7„,?;„,(J„))ieiK 

contained in A en and passing through (u„, j n , v n , S n ) where 

(u"(0),uf r (0),u?(0),^ r (0)) = (u n , ln ,v n ,6 n ) 

(cf., e.g., gDl Proposition 2.39]). 

In view of the regularity A £n C Pi (see Q3.4p ). we obtain the uniform bounds: 

e n \\u\W)\\ 2 + \\<(t)\\l + IK(*)IIW) + KWII2 + IKWII W) < O, (3.92) 
where the constant C > is independent of £ and e n . Now, for all T > 0, the functions u e ™ , ufp , u e t n , w^p 
and y/e^u% are, respectively, bounded in L°°(-T, T; -ff 2 (fi)), '(-T,T; ff 3/2 (r)), L°°{-T, T; H 1 ^)), 
L co {-T,T;H 1 / 2 {Y)) and L°°(-T, T; L 2 (ft)). Thus, there is a function u and a subsequence (not rela- 
belled), in which, 

m £ " ^ u in L°°(-T,T;iZ 2 (Q)) (weakly*), (3.93) 
u[f l ->> u ]r in L°°(-T, T; H 3/2 (T)) (weakly*), (3.94) 

u £ t n u t in L°°(-T, T; H 1 ^)) (weakly*), (3.95) 
uff r - ut|r in L°°(—T, T; i/ 1 / 2 (r)) (weakly*), (3.96) 

E n uf t " -> in L°°(-T,T; L 2 (ft)) (strongly). (3.97) 
The above convergence properties yield 

m £ " -> u in C(-T, T; ff^fi)) (strongly) (3.98) 
owing to the following embedding 

{u e £°°(-T, T; tf 2 (ft)) : w t e L°°(-T, T; ff 1 ^))} °-> C(-T, T; H 2 ^^)), (3.99) 

which is compact for every r\ £ (0, 1) (see, e.g., [37]). The strong property (|3.98[) allows us to identify 
the correct limit in the nonlinear term when e n — > 0. Moreover, from (|3.93[) and the fact that H 2 (f2) <^ 
C° (H), it follows that 

sup \\f{u^)-f(u)\\ 2 < sup Q4K"WL>WUIK"M-«WII 2 (3.ioo) 

te[-T,T] te[-T,T] 

<C(fi) sup ||u e "(t)-w(t)|| 2 , 
te[-T,T] 

for some positive (increasing) function : M + x M + — > R + , independent of n and e n . By virtue of p. 981) 
it is then easy to see that 

f(u^) -> /(u) in C(-T,T-L 2 (Q)) (strongly). 

It follows that u is a weak solution of the limit parabolic problem on M. In particular, (u n ,7n) = 
(w"(0),wp r (0)) -> (u(0),U|r(0)) in V 1 . Hence, we have that (u(0),«| r (0)) = («*,7*) 5 and therefore 
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(u(0),U| r (0)) G V 2 . As (u,u\ T ) is a complete orbit through (u*,t*), it follows that (it*, 7*) G Aq. It 

remains to show that v* = Au* — f(u*) and 5* = — <9„7* — 7*, in which case (it*, 7*, u*, <5*) G ^4o- 
Now by ([3U7|l and it follows that 

||e„u£(0)|| = V^\\V^< t (0)\\ < v^c, 

and so e n uf t (0) in L 2 (£7) as e„ —> 0. With this at hand, 

u?(0) = -e„u£(0) + Au"(0) - /K(0)) 

= -e B u£(0) + Au* - /(«*), 

so that 

<(0) Au* - f(u*) in L 2 (ft) (weakly). (3.101) 
Since u?(Q) = v n , then with (|3.101l) we have that 

v* = Au* - /(it*). (3.102) 

Similarly, since 

u^ r (0) = -^«"(0)-«f r (0) 
and since u" r (0) = 7* and it^ r (0) = S* , then 

5* = -3 n7 * - 7*. (3.103) 

We know (it*, 7*) G ^4o, so p.!02p and (|3.103j) imply that («*, 7*, v* , 5*) G i , in contradiction to IpTM]) . 
This proves the assertion and completes the proof. □ 

4. Exponential attractors 

Exponential attractors (sometimes called, inertial sets) are positively invariant sets possessing finite 
fractal dimension that attract bounded subsets of the phase space exponentially fast. It can readily be 
seen that when both a global attractor A and an exponential attractor M. exist, then A C A4, and so the 
global attractor is also finite dimensional. The existence of an exponential attractor depends on certain 
properties of the semigroup; namely, the smoothing property for the difference of any two trajectories 
and the existence of a more regular bounded absorbing set in the phase space (see, e.g., [20], [21]'). 

The main result of this section is the following. 

Theorem 4.1. For each e G (0,1], the dynamical system (S^,?^) associated with il. 1\) - [T73\) admits an 
exponential attractor M. e compact inH. e , and bounded inC £ . Moreover, there hold: 

(i) For each t > 0, S e (t)M e Q M e . 

(ii) The fractal dimension of M £ with respect to the metric T-L £ is finite, namely, 

dim F (Ms, H e ) < C £ < 00, 

for some positive constant C £ which depends on e. 

(Hi) There exist g > and a positive nondecreasing function Q £ such that, for all t > 0, 

<tist ne (S £ (t)B,M e ) < Q £ (\\B\\ He )e- e \ 

for every nonempty bounded subset B ofH £ . 

Remark 4.2. Above, 

a- 1 ka v \ v lnM« e (M E ,r) 

diniF(A4 e , He) ■= hmsup < 00, 

r^o - In r 

where, [iu s {Z,r) denotes the minimum number of r-balls from T-L £ required to cover Z C H £ . 
Corollary 4.3. There holds 

dim F (A £ ,H £ ) < dim F (M £ ,H £ ). 
As a consequence, A £ has finite fractal dimension which depends on e > 0. 

Remark 4.4. Unfortunately, we cannot show that the fractal dimension of Ai £ is uniform with respect to 
e > (see the subsequent lemmas). 

The proof of Theorem 14 . 1 1 follows from the application of an abstract result tailored specifically to our 
needs (see, e.g., [2TJ Proposition 1], [22], [30]; cf. also Remark [4.101 below) . 
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Proposition 4.5. Let (S £ ,T-L £ ) be a dynamical system for each e > 0. Assume the following hypotheses 
hold: 

(CI) There exists a bounded absorbing set B\ C T> £ which is positively invariant for S £ (t) . More 
precisely, there exists a time ti > 0, which depends on e > 0, such that 

S e {t)B\ C B\ 

for all t>ti where B\ is endowed with the topology ofTL £ . 
(C2) There is t* > t\ such that the map S £ (t*) admits the decomposition, for each e £ (0,1] and for 
all ipo, 6*o <G Bl, 

S e (t*)vo - S £ (t*)9 a = L e {<p , B ) + Re(<Po,6 ) 
where, for some constants a* <S (0, |) and A* = A*(il,t*) > with A* depending on e > 0, the 
following hold: 

\\L e (ip o ,0o)\\n,<a*y o -e o \\ He (4.1) 

and 

\\Re(^o,0 Q )\\ Ve < A*\\cp Q -0 o \\ Hc . (4.2) 

(C3) The map 

(t, U) i — ^ S £ {t)U : [t*,2t*] x Bl -> Bl 
is Lipschitz continuous on Bl in the topology of TL £ . 
Then, (S £ ,H £ ) possesses an exponential attr actor M. £ in B x £ . 

We now show that the assumptions (C1)-(C3) hold for (S £ (t) , T-L £ ). We begin with a higher-order 
dissipative estimate in the norm of T> £ . 

Lemma 4.6. Condition (CI) holds for fixed e € (0, 1]. 

Proof. The proof is very similar to the proof of Lemma [3.161 Indeed, let e € (0, 1], ipo = (ito,wi) € T> £ 
and f (t) = S £ (t) (po- In this setting, we differentiate (|1.1[) - (|1.3[) with respect to t and let h = u t - We set 
(3 in (|3.28p to be f3 = $ where we recall that d > is due to assumption (| 1 . T|) . Then we easily obtain 
the analogue of the differential inequality f|3.6T[) except that the size of the initial data now depends on 
the norm of T> £ , i.e., tpo = (uq,ui) € T> £ (here the initial conditions are not necessarily equal to zero). 
Thus, after applying (|3.60[) and (|3.61[) . there exist a positive and nondecreasing function Q and a constant 
C > such that 

||(M*), ht(t))\\n c < 0(11 (M0), fh(0)) IkJe"^*/ 2 + C (R) (4.3) 
(Q, loq and C arc independent of e) with R > such that H^oll-He < R- Arguing as in Theorem 13.161 by 
exploiting i? 2 -elliptic regularity theory, we also deduce 

||V (t)fv e < Osd^olbJe-^*/ 2 +C(R), (4.4) 
for some new function Q £ which depends on e > 0. Indeed, using the equations (|l.ip - (|1.3[) . it is not 
difficult to show that there holds || (h(0), ht(0)) \\u s < "^llvollc whence (|4.4[) . Consequently, there 
exists R± > (independent of time, e > and initial data) such that S £ (t) possesses an absorbing ball 
Bl = P>T> e (Ri) of radius R\ centered at 0, which is bounded in T> £ . This establishes condition (CI). □ 

Remark 4.7. Unfortunately, the bound in the space T> £ is not uniform as e — > + . Indeed the function 
Qe (•) in (|4.4p blows up as e — > + . Finally, arguing in a standard way as in Theorem 13. 171 Bl is in fact 
exponentially attracting in T-L £ . 

Lemma 4.8. Condition (C2) holds for each fixed e £ (0, 1]. 

Proof. Let e € (0,1]. Let (fio,0o £ B\. Define the pair of trajectories, for t > 0, ip(t) = S £ (t)ipo = 
(u(t),u t (t)) and 0{t) = S £ (t)6 = (v(t),v t (t)). For each t > 0, decompose the difference ({t) := ip(t)-9(t) 
with Co := tfo — 9q as follows: 

where (p(t) = (u(t),iit(t)) and 0(t) = (v(t),v t (t)) are solutions of the problems: 

eu t t + u t — Aw = in (0, oo) x 

d n u + u + u t =Q on (0, oo) x T (4.5) 
<p(0) = tfio — in fl 
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and 

ev tt + v t - Av = f(v) - f(u) in (0, oo) x fl 

d n v + v + v t = on (0, oo) x T (4.6) 

0(0) = in fi. 

By estimating along the usual lines, multiplying (|4.5jl i by 2u t + u in L 2 (f2), we easily obtain the 
differential inequality, for almost all t > 0, 

+ uj 7 N € < 0, (4.7) 

for some positive constant uij sufficiently small and independent of s, and for 

N e =N e {(p{t)) :=e|l^(*)l| 2 +^W^(*)) + l|Vu(t)|| 2 + ||u(t)||| 2(r) . (4.8) 

Obviously, J\f e is the square of an equivalent norm on H e , i.e., there is a constant C > 0, independent of 
£, such that 

C- l \W\\u^<K^)<C\m 2 He . (4.9) 
Following (|4~7f and (|4.9j) . we have that, for all t > 0, 

11^)11^ <C||^||^e-^. (4.10) 

Set i* := max{£i, ^ In (4(7)}. Then, for all t > t* , (JUTJ holds with L e = (p(t*) and 

a* = Ce-" 7t * < -. 

2 

We now show (|4.2[) holds for some A* > 0. First we observe that 

2(f(v) - /(«), w«) £a( r) = ^2</(w) - /(«), « t ) L 2 (r) - u 

- 2<(/» - /'(u))w t) C t ) i2(r) + 2{f(u)z u v t ) L 2 (T) . 

Next we differentiate the second equation of (|4.6[) with respect to i, multiply the first equation of (|4.6[) 
by 2(— A)u t in L 2 (f2) and insert (|4.11[) into the result to produce the differential identity, which holds for 
almost alH > 0, 

^ {elMli + INI W) + II Av\\ 2 + 2(f(u) - f(v),Vth'(T)} + 

+ 2£||^ t |li2 (r) +2||v t || 2 (412) 
= 2((/» - /'(u))V«, V«t) - 2(/'(u)Vz, V«*)+ 

+ 2(/(«) - f(u),v t ) L 2 (r) - 2((/» - f{u))v u v t ) LHr) + 2</ / (u)*, v t ) L , [T) . 

Recall that z :— u — v denotes the difference of any two weak solutions of (|1.1[) - (|1.3[) and is estimated in 
(|3.8[) . Arguing, for instance, as in [23l (6.11)-(6.13)], we estimate the products on the right hand side of 
dHHl), for all * e (0,**), using (fTTlj) . Lemma ESI and the embedding H 2 (fl) ^ C°(U), as follows: 

2|<(/» ~ f'(v))Vv,Vv t )\ < C{\ + \\u\U + \\v\U) WzUvUVvtW 

<c E (t*)\Mk + \\\Vvt\\ 2 , 

2\(f'(u)Vz,Vv t )\ <C(l + !| M || 2 )||Vz||!|V« t || 

1 (4-14) 
<C e (**)HCo||| £e + i l|Vi? t || 2 , 



2\(f(u)-f(v),v t ) L 2 { T)\<C\\z\\i\\vt\\L* { r) 

<cn( \\ 2 H E + \\\vt\\h {r) , 

2 I((/'( U ) - f'(v))v t ,vt)LZ(T)\ < (7(1+ ||u||c°(r) + Nlc°(r)) \\z\\c«{r)\\vt\\mr)\\vt\\mr) 



(4.15) 



<^*)llCo||^ + i ll^||i 2( r ) 



!„_„, ( 4 - 16 ) 
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and 

2 \(f'( u ) z t,vt}L2(r)\ <C(l + ||w||co(r)) IMU 2 (r) INU^r) 

1 I (4.17) 
<C(r)||Co|lL + 4ll^H^(r)- 

We emphasize again that by Lemma 14.61 the constants C = C £ (t* ) in estimates (|4.13|) and (|4.14|) depend 
on e > 0. After combining (|4.13[) - (|4.17[) with the identity (|4.12p . we arrive at the differential inequality, 

^ {^IMIi + INI W) + IIA^II 2 + 2{f(u) - f(v),v t ) LHr) ] < C E (n\\Co\\ 2 n e (4.18) 

(Recall that by the definition of a strong solution in Definition 13.51 v u G L 2 (0, oo; L 2 (T))). Now by 
integrating (|4.18|) over (0,t*) and once again applying the estimate (|4.15|) . we are left with the bound 

^(t*)lli + l|A«(t*)|| 2 <C e (i*)llCo||^. 
By standard -ff 2 -elliptic regularity estimates (see (|3.70[) and (|3.7ip above), we obtain 

\mn\\v e <c s (tnuo\\n e - (4.19) 

Inequality (|4.2j) now follows with R £ = 8(t*) and A* = C £ (t*) > 0. This finishes the proof. □ 
Lemma 4.9. Condition (C3) holds. 

Proof. We proceed exactly as in the proof of Lemma l4~6l differentiating (|1.1|) - (|1.3|) with respect to t and 
letting h = Ut- This time we obtain the bound 

\\vt(t)\\u e < Qe{R) 

for (ft = (u t ,Uu), and some function Q £ , depending on e > 0, where the size of the initial data now 
depends on the norm of B\. Hence, on the compact interval [t*,2i*], the map t i— > S e (t)(po is Lipschitz 
continuous for each fixed (po £ B\\ i.e., there is a constant L £ = L e (t*) > (which depends on e > 0) 
such that 

\\S e (h)v - S e {t 2 )w\\u c < L e (t*)\h - t 2 \. 
Together with the continuous dependence estimate p.8j) . (C3) follows. □ 

Remark 4.10. According to Proposition l4.5l the scmiflow S £ : H £ —> H £ possesses an exponential attractor 
A4 £ C Bl, which attracts bounded subsets of B\ exponentially fast (in the topology of H £ ). In order 
to show that the attraction property in Theorem 14.11 (iii) also holds, we can appeal once more to the 
transitivity of the exponential attraction [22l Theorem 5.1] and the result of Theorem 13.171 (also see 
Remark El} . 

In contrast to the standard case of Dirichlet boundary conditions, where we have a complete treatment, 
due to References [22] [41], the situation with boundary condition (|1.2[) remains essentially less clear. The 
following important questions remain open: 

• Higher-order dissipative estimates which are uniform with respect to e > 0. 

• Finite-dimensionality of the exponential attractor M £ (and global attractor A £ ) which is uniform 
in e > 0. 

• Existence of a robust (Holder continuous in e € [0, 1]) family of exponential attractors {.M £ } . 

5. Appendix 

To make the paper reasonably self-contained, we include the statement of a frequently used Gronwall- 
type inequality [441 Lemma 5]. 

Proposition 5.1. Let A : R + — > M + be an absolutely continuous function satisfying 

— A(t) + 2r?A(t) < h(t)A(t) + k, 
dt 

where r\ > 0, k > and J* /i(r)dr < rj(t — s) + m, for all t > s > and some m > 0. Then, for all t > 0, 

kp m 

A(i) < A(0)e m e-"* + — . 
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